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Study guide 


Before reading this unit you should read those sections of Unit 25 which refer to 
the calculation of partial derivatives. 


You should read Sections 1-3 before watching the television programme, but you 
can choose to read Section 4 either before or after the programme. 


1 Introduction and definitions 


1.1 Introduction 


This unit is devoted almost entirely to what can truly be described as one of the 
“gems’ of late eighteenth and early nineteenth century mathematics—the solution 
of partial differential equations with given boundary conditions—although in the 
space available we shall be able to examine only one simple case in any detail. 


My objective will be to give you a glimpse of the power of the methods which you 
have studied so far —and incidentally, and more prosaically, to enable you to 
revise some of the topics of this course in what I hope is an interesting fashion, 


The subject which we shall study arose in a period when a number of great 
mathematicians were turning their minds to the investigation of various questions 
arising in what we should call applied mathematics. To their credit, the men of this 
age would have drawn no clear distinction between what was ‘pure’ and what 
‘applied’; they were simply intent on obtaining answers to interesting questions, 
and in searching for these answers they incidentally laid the foundation for 
modern mathematical analysis as well as for modern physics. 


A catalogue of the names associated with the subject reads like an honours list: 
Euler, Laplace, Lagrange, Poisson, Cauchy, the Bernoullis, Fourier, and many 
others, Before we become involved in the details of the single example which we 
have time to study, it will perhaps give you a deeper understanding of the 
relevance of the topic if I say just a little about the range of subjects which these 
men considered. 


A number of apparently distinct areas of. physics — electrostatics, magnetism, 
hydrodynamics, aerodynamics, heat conduction and wave motion— use models 
which involve differential equations with partial derivatives. Such an equation is 
called a partial differential equation. The equation 

#0 00 

ax?" at 
is an example of a partial differential equation, and it can be used to model the 
flow of heat along a thin metal bar; in this case 0(x,1) is the temperature at 
Position x at time 1. Having deduced that @ satisfies such an equation, we are very 
little further forward unless we are able to determine 6 explicitly, using any 
available information about the temperature at some particular point in time and 
how heat may be lost from the system. 


The break-through in such problems came with the development of trigonometric 
series and the publication of Joseph Fourier’s Theorie Analytique de la Chaleur 
(1822), which is the source of many important methods in mathematical physics, in 
particular for problems involving partial differential equations under given 
boundary conditions. Certainly the subject had been advanced by Euler, 
D’Alembert, David Bernoulli and others, but it was Fourier who made the 
situation absolutely clear by establishing the fact that an ‘arbitrary’ function (in 
fact a piecewise continuous function) could be expressed in the form 

D (A, cos nax + B, sin nax). 

n=0 
Despite the support of Euler and Bernoulli, the idea was so original when first 
presented in 1807 that it raised some vigorous opposition, notably, it is said, from 
Lagrange (the leading French mathematician of the day). The reluctance of some 


See, for example, Chapter 22 
of Mathematical Thought from 
Ancient to Modern Times, by 
Morris Kline (Oxford 
University Press, 1972). 
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mathematicians to accept Fourier’s results must certainly have been due in part to 
the fact that it called into question the very meaning of the fundamental concept 
of a function. Ultimately, the ‘Fourier series’ was universally adopted, and the 
questions which it raised spurred mathematicians to re-examine the foundations of 
the calculus and to restore it in its modern form as mathematical analysis 


The essential elements of Fourier’s method are easy to describe. First derive the 
appropriate partial differential equation, then assume that the required function 
can be expressed in terms of a Fourier series, and finally use the equation and the 
initial conditions to determine the constants in the series. This is precisely the 
outline of the core of this unit; but, before we start on this programme, I should 
like to look at the problem that we shall attempt to solve in more detail. 


Imagine a coiled spring stretched tightly between two fixed points A and B 
(Figure 1). Now suppose that a point in the middle of the spring is pulled to one 
side, towards B, and then released, while the ends are kept fixed (Figure 2). The 
spring will start to vibrate, with each point of the spring moving to and fro in 
some fashion. Our task is to determine how each point will move, and to discover 
where it will be at each subsequent instant of time. In our solution of this problem 
there will be three main factors: 


(i) the physical properties of the spring, which will give rise to an equation 
known as the wave equation; 


(ii) the physical constraints applied to the spring, namely that it is fixed at 
both ends, which give rise to the so-called boundary conditions; 


(iii) the initial behaviour of the spring, in other words, what we do to the 
spring in order to make it move. This initial behaviour produces the initial 
conditions for the problem. 


The exact meaning of these terms will become clear as we proceed through the 
unit. 


Although I have chosen to focus on the motion of a vibrating spring, the methods 
which we shall investigate are very general, and apply to the equations used for a 
wide variety of physical situations including heat flow, the sound produced by 
musical instruments and the motion of fluids. 


The essence of each of these physical situations is that it can be modelled by a 
partial differential equation, and so | must explain in more detail what | mean by 
such an equation, 


1.2 Partial differential equations 
Previously in this course we have solved a differential equation such as 


a " dy 
oa + oe +ey=0, 
which is an example of an ordinary differential equation, and is characterized as 
such by the fact that the unknown function y(x) is a function of the single variable 
x. We are now intending to discuss differential equations which apply to functions 
of two variables, Such a function, U(x,f) say, can vary both because of changes in 
x and because of changes in 1. 


a 


As an illustration, consider the equation oe = . The function U(x, t) = e*' is 


Ox 6t 
a solution of this equation, as we can easily verify. We have 
eu Luseees 
ae te“ and aoe xe", 
and it follows that xf = na and so 
Ox et 


U(x, =e 


AAW - BB 


Figure 1. A coiled, stretched 
spring. 


ASIANA BB 
— 
Figure 2. The same spring 


with the middle pulled 
towards B. 


We shall see later how this use 
of the term ‘boundary 
conditions’ relates to the use 
in Unit 6, Section 3. 
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is a solution of the partial differential equation 


eu _ aU 
oa Fie 
As you may have guessed, the general problem which we should like to solve is to 
discover all of the functions which satisfy a given partial differential equation, and 
this is the process which we call solving the equation. Generally speaking, solving a 
partial differential equation is by no means easy, and the following illustration 
may help you to understand why. 


Suppose that we are given the rather simple looking partial differential equation 
OU(xt) _ 1 
oa ta 


The left-hand side is the partial derivative of U(x,t) with respect to x keeping 1 
constant, so to find U we may keep 1 constant and integrate with respect to x. We 
then obtain 


U(x, t)h=x+C. 
However, the constant C is only constant if r is kept fixed, and for different values 
of ¢ the value of C will be different, so that C is a function of 1, Thus the most 
general solution is 

Ulex) =x +f) 
where f is an arbitrary function, You can verify that any such function U(x, 1) is a 
solution by differentiating partially with respect to x. 
So you seé that while the solution of an ordinary differential equation will involve 


arbitrary constants, the solution of a partial differential equation will involve 
arbitrary functions. 


As we increase the order of the partial derivatives in a partial differential equation, 
so we introduce more arbitrary functions into the solution. Thus the partial 
differential equation 

a2. U 

GU Lo 

Oxét 
might be expected to contain two arbitrary functions in its solution, and the 
following argument shows that this is indeed so. 


Integrating with respect to x (for constant 1) we see that > = f(t), where f is an 


arbitrary function (of one variable). Now integrating this equation with respect to 
t (for constant x), we get 


U(x.) = frou + a(x), 


where g is another arbitrary function. This solution can be written in the form 
U(x, 0) = h(t) + g(x) 
where both h and g are arbitrary functions. 


In order to make any progress it will be necessary for me to restrict the class of 
partial differential equations to be discussed. Partial differential equations are 
classified in much the same way as ordinary differential equations (see Unit 6), so 
that a linear second-order partial differential equation in two independent variables 
will be one of the form 


U @u 6U eu 
at (x,t) a + d(x,t) a e(x,0) = + f(x. NU = g(x,t) 


@U 
6x? 


Pd 
a(x, 0) + B(x, 1) = 


where a, b, c, d, e, f and g are known functions of x and 1, and at least one of a, b 
and c is non-zero, 


tno longer 
fixed 
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This is still far too wide a class to tackle in one unit, and in fact our concern will 
be with equations which (in addition to obeying the above conditions) are 
homogeneous (i.e. g(x, f) is the zero function) and constant-coefficient (i.e. a, b, c, d, 
e and f are all constants). Thus, our attention will be confined to equations of the 
form 

eu @U @U au. au 

Be ae ae ag ae TIA a) 


where a, b, c, d, ¢ and f are constants. 


a 


These are the analogues of the ordinary differential equations of the form 
dy 
dx? 
which are discussed in Section 1 of Unit 6. Theorem | of that section also has the 


following analogue for partial differential equations, which is of extreme 
importance to us for our work later in this unit. 


Theorem 1 


If wand v are two solutions of a linear homogeneous partial differential equation, 
then for all values of A and B, the function Au + Bo is also a solution of the 
equation. 


dy 
a +b +cy=0 


Proof (for the second-order constant-coefficient case) 
Let U = Au + Bo, where u and v are solutions of Equation (1), Then 


@U e 
CN = ana lau + Bo) 
Ou &v 
=aA ie + Boa 


@u ev 
=A (053) 2 8 (053) 
and similarly for the other terms in the equation. Thus 
@uU &@uU @uU ou ou 


get aa eee ae TE 


a 
Ou ry 
=A (oS + +h) +0 (055 + +) 
=0 (since u and v are solutions of Equation (1)), 
so that U is a solution of Equation (1). 


In fact this theorem holds (as stated) for all linear homogeneous partial differential 
equations; the general proof is longer than the above but very similar in its 
method. 


Before we go any further, I must make sure that you are confident with partial 
derivatives, so if you have any doubts about your ability try the following 
exercises. 


Exercise 1 


5 git: oF : ‘ 
Find ae and = for each of the following functions F(x, y} 


() Fes) =* 

y 
(ii) F(x.y) = Ax? + 2Bxy + Cy? 
(iii) F(x, y) = sin(2x + y) 


(iv) F(x, y) = @(2x + y), where ¢ is an arbitrary differentiable function of one variable. (Use 
the notation ' for the derived function of @.) 


(v) (x,y) = (ax? + by?) where ¢ is an arbitrary differentiable function of one variable. 
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Exercise 2 
Write down simple expressions for the following: 


tim L(x + 6x) — f(x) 


x0 ox 
(i) lim F(x + dx,1) — F(x,t) 
x0 ox 
a 
Lear + 6x,1) — oF oe of) 
ati 
x0 ox 7 
Exercise 3 


OF OF PF OF é . 
ae Dyax" axdy a for each of the following functions F(x, y). 


(i) Fup) = xy 


Find 


(i) F(xyp) = 38 + 2x7p + y? 
(iii) Foxy) =e*” 
Exercise 4 


Verify that each of the following functions is a solution of the corresponding partial 
differential equation. (f and g are arbitrary differentiable functions of one variable, and a 
and b are given constants.) 
; @u _eueu 
(i) U=fixgu); pes la Op 
aU . aU 


(ii) U = fay — bx); ois yr Se 


(iii) y = f(x — 0) + glx + 0; 


Exercise 5 
Determine which of the following partial differential equations are homogeneous linear, and 
which are constant-coefficient. 
OF OF 
Pe aera ss 
(i) x me +P SR 0 


(ii) Fees ee ee 
Gt + * ae oy THF 
ar ey) oF oF 


FOF LaF = oy 
ax? * ay? S 


(iii) x2 | aa 


.) OF OF . 
(iv) Pe ied (where a is a real constant) 


[Solution to Exercises 1-5 on p. 35] 


Summary of Section 1 


1. In this section you have seen examples of partial differential equations and their 
solutions, and you have seen that such solutions involve arbitrary functions, as 
opposed to the arbitrary constants which occur in the solution of ordinary 
differential equations. 


2, In this unit we shall restrict our attention to homogeneous, constant-coefficient, 
linear, second-order partial differential equations; ie. equations of the form 

@U #U @U 6U 6U 

Be eget oe: ea eg TIO 


where a, b, c, d, e and f are constants. 


a 


3. Theorem 1: If u and v are two solutions of any homogeneous, linear partial 
differential equation, then for any constants A and B, the function Au + Bo is 
also a solution of the equation, 
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2 Properties of coiled springs and elastic 
strings 
2.1 Hooke’s law 


A coiled spring is a device, usually made of steel, which can be altered in length, 
but which returns to its original shape when any external forces are removed. Such 
a spring may be extended, in which case we say that it is in a state of tension, or 
possibly compressed, which we call a state of compression. Our objective in this 
section is to model the behaviour of such springs in some detail, and in order to 
do this we shall need to make certain simplifying assumptions. 


The springs which you have met previously in this course were assumed to be ideal 
springs, that is to say, they were assumed to obey Hooke’s law and to be light. (By 
‘light’ we mean that the mass of the spring can be ignored because it is small 
relative to the other quantities under discussion.) You will recall that Hooke’s law 
states that an ideal spring when uniformly extended will exert a force F, given by 


F=K(I—Ip) 
where K is the stiffness of the spring, /o is its natural length and | is its extended 
length. We shall find it convenient to replace K by ‘t and the constant k is then 
0 


known as the modulus of elasticity of the spring, or more briefly its modulus. We 
then have 


Peony 
lg 


The springs under consideration in this section are no longer assumed to be light, 
but I shall make the assumption that they continue to obey Hooke’s law. (This is 
an assumption which can be validated only by experiment.) In future in this unit 
all references to ‘springs’ are intended to mean heavy springs which obey Hooke's 
law. 


Later in this section I propose to model the behaviour of a spring which is 
originally stretched horizontally between two fixed points. The spring is then 
disturbed slightly from its rest position, and our task will be to predict its 
subsequent motion. The essential point to note at this stage is that our starting 
point could well be a spring which has already been extended, but before I 
continue with this line of thought I must define the term tension in the spring. 


We consider a stretched spring with ends A and B which may, or may not, be 

stationary, and consider an arbitrary point P along its length (see Figure 1). —T 
The section of spring to the right of P exerts a force towards B, of magnitude T 
say, on the left-hand section; and from Newton's third law we deduce that the left- A-WAWANAAANA 8B 
hand section applies an equal and opposite force on the right-hand section. The T— 

magnitude of this force is called the tension in the spring at P. The value of T will Figure 1. The forces acting at 
in general depend upon the choice of the point P, but we shall assume throughout P on the left-hand and right- 
our discussion that the spring is always stretched rather than compressed, so that hand sections of a stretched 
T is always positive. oe 


My first task is to examine the relationship between the tension in the spring and 
its physical shape. 


The springs which I shall consider will be long compared to their width, so that in 
future diagrams I shall represent them by straight lines. 


Suppose that the spring is now stretched between two fixed points A and Ba 
distance L apart, and that the spring is stationary. This situation will in future be 
known as the reference state of the spring. 


While in its reference state the tension at all points of the spring is the same, Ty Ta Th 
say, for each small element of the spring must experience equal and opposite forces 
at its ends, and therefore the tension must be the same at both ends of the Pow 

element. (If this were not so, the element would begin to move in the direction of Figure 2. The forces at P and 
the greater force.) Q must be equal. 
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We now choose an arbitrary point P on the spring (and mark it say with a white 
dot), and our objective will be to describe how P moves when the spring is 
disturbed. Suppose initially that while the spring is in its reference state the point 
P is at a distance x from A. Suppose also that Q is a point close to P, at a 
distance x + 6x from A while the spring is in its reference state. 


The element of the spring between P and Q of length 6x had an original 
unstretched length of 5a, say, so that applying Hooke’s law to this small section of 
spring we obtain 
[ix — da 
Py 


We can solve this equation for da and thus find an expression for the unstretched 
length of the element PQ in the form 


a= { Jo (l) 


k 
To + k 
Let us now consider what happens if the spring is set in motion by a small 
disturbance along its length. Suppose that at time 1 the point P on the spring has 
moved so that its distance from A is s, and further suppose that Q is a distance 
s + 6s from A. Using Hooke’s law again we have 


Tx (5%), Q) 


since the length of the element PQ is now és, and T is the tension at P. (This 
equation is approximate since we assume that the spring is uniformly extended 
between P and Q in order to apply Hooke’s law.) 


Now, eliminating da between Equations (1) and (2), we obtain 
6s 
T=(To+ Li -k, 


so that in the limit as 6x + 0 we eliminate the local effects of non-uniform 
stretching of the spring and obtain an exact result 


T=(h+ he k. (3) 
Ix 


d: 
(For the time being I am regarding ¢ as fixed, so that it is unnecessary to use a 
partial derivative in Equation (3).) 


In Equation (3) you should notice that “ is a quantity which depends upon the 


way in which the spring is distorted, so that this equation relates the geometric 
shape of the spring to the tension at any point of the spring, Also I have not 
actually used the condition that A and B are fixed, so that Equation (3) would 
also apply in the case where A and B are free to move. The above equation is used 
in the following example, which is intended to get you accustomed to the idea of 
tension varying along the length of a spring. 


Example 1 
A stationary uniform spring of mass Mp, modulus k and unstretched length / 
hangs vertically and supports a weight of mass M. Find the length of the stretched 
spring. 
Solution 
We are told that the spring is uniform, which indicates that its mass per unit length 
M 
when in its unstretched state is constant and equal to = (When stretched, 
0 
however, the mass per unit length is no longer constant because the stretching is 


not performed evenly along its length.) The solution now falls into two parts: an 
application of Newton’s second law, then an application of Hooke’s law. 


reference state 


A PO 

1 ede Hl 

\ bo | 
1 \ 

A *  P&sQ B 

——L. — 


disturbed state 


Figure 3. The reference state 
and a disturbed state of a 
stretched spring. 
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(i) Apply Newton's second law to a small element of the spring 


We choose the variables as shown in Figure 4, and consider the forces acting on a 
small element és of the stretched spring. 


The left-hand side of the figure shows the spring in its reference state (with gravity 
‘switched off’) which corresponds to the case in which Ty = 0 and the spring is 
unstretched. This element was originally of length 5x, and its mass is therefore 
Mo 

is: 6x. 


Suppose now that the tension in the stretched spring is Tp at the point P and Tg 
at the point Q. The element PQ experiences three forces (see Figure 5): 


(i) its weight downwards, 
(ii) the tension at P upwards, 
(iii) the tension at Q downwards, 


giving a net downwards force of 


Fert ey. 
Io 


Since we are assuming that the stretched spring is stationary, the acceleration of 
the element PQ is zero, and from Newton's second law we have 
force = mass x acceleration, 


so that 
To — Te + MS 5x a0. 
lo 
We now write the difference in tension between Q and P as 67, ie. 


6T = To — Tp 
It follows that 


ook 
eet a 
and in the limit, 
dT Mog 
os ae: 4) 
dk lo @ 
Now integrating with respect to x, 
Mog? 
Te) = -—" +6, (5) 
Io 


where C is a constant. T(x) is thus the tension at the point P which was at a 
distance x from A when the spring was in its reference state. 


In order to determine C we consider the forces acting on the mass M: 


(i) a force upwards of Ty = T(Io) = — Motte +C=C— Mog, 

(ii) a force downwards of Mg; 

so that the net force downwards is 
Mg + Mog —C, 

and this is equal to zero since the mass M is stationary. It follows that 
C=(M+ Mo). 


(You should notice in passing that we could also have found the value of C by 
examining the tension at the point A and equating this to the total weight hanging 
from A.) 


Here, unlike our previous 
situation, the “unstretched 
state’ is our ‘reference state’. 


unstretched > 


spring ~ 
(with gravity B 
‘switched ol") etched 
spring 
(with gravity 
“switched on’) 


Figure 4. The spring in its 
unstretched and stretched 
states. 


Mis ay 
a Te 


Figure 5, The forces acting on 
PO. . 


Ty (= Ths) 


Me 


Figure 6. The forces acting 
on the mass M. 
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Equation (5) now becomes 
T(x) = Moet +(M + Mo)g. (6) 
lo 


Thus Newton's second law tells us how the tension T(x) depends on x. 
(ii) Apply Hooke's law 


From Equation (3) on page 10 we have, since Ty = 0 in the reference state we have 
chosen, 


ds 
T(x) =k (@ -1). (7) 


Our problem is to find the length of the stretched spring, which is s(/o) = / say, 
and this would be solved if we could find the function s(x). Eliminating the 
unwanted variable T(x) between (6) and (7), we obtain 

ds Mogx 


pe tl) 


+(M + Mo)g. 


4 ds 
and solving for EY have 
a” Re k 
Integrating with respect to x, we obtain 


Mogx? _ [(M + Mo)g 
2klo + eee x+y 


for some constant C,. 


ds Mogx o (M + Mo)g v1. 


s=- 


But s = 0 when x = 0, and therefore C, = 0, so that 
2 
Besar (t+ Mow Mole i} (8) 


We ase asked to find s(/o), so substituting x = /p into (8) we have 


Mogi. [(M+M 
at) = — Hoe ; og 


= (te ion + hols + 1) 


te i)' 


When you derive an answer to a problem like this it is often useful to check its 
validity for extreme cases. Thus, for example, we see that if k is very large, 
corresponding to the case of a very stiff spring, then from Equation (8) we obtain 
s = x, which is just what we should expect. 


2.2 Vibrations of a stretched spring 


We now come to an important point in this unit, namely the point at which we 
define the problem which we shall eventually solve. 


We consider a heavy spring, at first in its reference state, stretched horizontally 
and tightly between two fixed points A and B at a distance L apart, and then 
disturbed by small longitudinal forces which impart a small motion to the spring. 
More specifically we make the following assumptions. 


(i) The motion of any section of the spring is small compared to the spring's 
length, takes place entirely along the length of the spring, and the end 
points are fixed. 


(ii) All frictional forces, such as that due to air resistance, may be neglected. 


(iii) The forces due to the tension in the spring act along its length and obey 
Hooke’s law. 
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(iv) The gravitational force on the spring is small relative to the tension, and 
so may be neglected. 


(v) The diameter of the spring is small compared to its length, and so may be 
neglected. 


Suppose that we choose an arbitrary point P of the spring, whose reference 
position is at a distance x from A, and suppose that the spring is disturbed slightly 
and begins to move. Our objective is to determine the position of the point P at 
any subsequent time. More specifically, suppose that the point P has moved to a 


position at a distance’s from A at time 1, then s = s(x, 1) is some function of x and reference state 

t, and our task will-be to determine this function. Suppose also that the point Q A PQ B 
close to P, whose reference position is at a distance x + 6x from A, has moved to 1 x lay 1 

a position at a distance s + ds from A at time 1 (Figure 7). 1 \ ‘ i 

Before the spring began to move it was in its reference state stretched uniformly A Ss P&Q B 
between A and B, and thus its mass per unit length was a constant, m say. (This is —L——— 

not equal to the unstretched mass per unit length.) disturbed state 


Figure 7. The spring in its 


In order to model the motion of the spring we need to bring together two Resice stain Geld caus titbed 


ingredients, as in Example 1: state. 
(0) Newton's second law (which specifies the acceleration if the force is 

known), 
(ii) Hooke’s law (which determines the tension). 


We let Tp be the tension at P and Tg be the tension at Q, at time 1. 
(i) Apply Newton's second law to a small element of the spring 


The small element of the spring between P and Q was originally of length 6x, so Tr Te 


that its mass is mdx. This element is acted on by two opposing forces, the tension Pt Q 
at P and the tension at Q. The centre of mass of this element is somewhere Figure 8 The forces acting on 
between P and Q and, since the element is small, its acceleration is approximately ‘a small element of the ares 


2 
a (x0) to the right. 


(This corresponds to the acceleration of a particular point of the spring, ic. a fixed 
value of x, and hence we keep x fixed in the differentiation.) From Newton's 
second law we have 


force = mass x acceleration, 


so that 
és 
To — Tr = mix 3, 
and therefore 
ae 
6x ar” 
In the limit as 6x + 0 the centre of mass of the element becomes P, so that 
Of, eae 
CLEP IC 9) 
ax "OF M2 


(ii) Apply Hooke’s law 
We now use Equation (3) on page 10 (derived from Hooke’s law): 


ds 
T=(hh +h) —k. 


The derivative in this equation refers to a fixed value of t. Now that we are 
c a SRS as 4 5 
treating ¢ as a variable we must write x to indicate that f is held constant during 


the differentiation, and we then have 


Ta(h+ =k (10) 
ox 
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Differentiating (10) partially with respect to x, we obtain 
(11) 
We wish to discover how s depends on x and t, and therefore we want to eliminate 


the superfluous variable T. We can do this by equating Equations (9) and (11) to 
obtain 


&s Cm 
a os 2 
maz = (To+ ky (12) 
k k 
The value of tks is positive, and it is common practice to put c? = bit to 


produce the partial differential equation 


&s 


(13) 


This equation is known as the (one-dimensional) wave equation, and is of crucial 
importance for the rest of this unit. In fact it is this partial differential equation 
that you will learn how to solve. The constant c is called the wave speed. 


The function s = (x,t) determines the distance of a point P from A at any time f; 
but it is more common to write the equation in terms of U = s — x, so that 
U = U(x,t) determines the distance of P from its starting point (see Figure 9). 


If we substitute s = U + x into Equation (13), we obtain 


(14) 


Exercise 1 

Verify that the substitution of s = U + x into Equation (13) gives Equation (14), 

[Solution on p. 35} 

The advantage of using U rather than s is that the conditions which follow from 

the assumption that the ends of the spring are fixed are quite simply stated as 
UO.) =0 
U(L,) =0 


These conditions are called the boundary conditions for the partial differential 
equation, and will be used extensively in Section 3. 


for all 1. 


2.3 Elastic strings 


An elastic string is just what it says it is, a string with elastic properties; but as a 
mathematical term it specifies that the string satisfies Hooke’s law. Springs and 
strings are therefore identical in behaviour except in one important respect; a 
string cannot be in a state of compression. 


We introduce the concept of an elastic string simply to give you another 
illustration of the use of the wave equation, and one which is perhaps a little easier 
to visualize. We suppose that the string is initially stretched tightly between two 
fixed points A and B, and then disturbed slightly in a direction perpendicular to 
the length of the string (see Figure 10). 


If P is a point on the string at a distance x to the right of A, then it can be shown 
that U(x,#), the perpendicular displacement of P from the line AB, will satisfy the 
wave equation 

cee 

Ox? ax e or? 


reference state 


y yams B 
| 1 

i a 

\ |g P. 

4 Us B 
—— Hs 


disturbed state 


Figure 9. Derivation of the 
function U(x,1). 


Figure 10, An elastic string 
vibrating. 
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where c is some constant, dependent on the physical properties of the string and 
its initial tension. (I ask you to accept this on trust since | feel that it would 
take too much of your time to give yet another derivation of the wave equation.) 


In this case the boundary conditions are again 
U(0,t)=0 
U(L,t)=0 


As in the case of a spring, the task is to determine the function U(x,r) from the 
initial behaviour of the string. 


hoe all ¢. 


In order to distinguish the two forms of motion, vibrations along the length of the 
string or spring are called longitudinal vibrations and vibrations across a string at 
right angles to its length are called transverse vibrations. 


Summary of Section 2 


1. In this section we modelled the behaviour of a ‘heavy’ spring, and we used 
Hooke’s law to derive the equation 


ds 
T=(Ih+ KE 
2. By combining this result with Newton's second law we showed that with certain 
assumptions the longitudinal displacement U of a stretched spring satisfied the 
wave equation 


for some constant c, where U satisfies the boundary conditions 
U0, =0 
viLD= of er allt. 


3. We introduced the terms longitudinal and transverse vibrations, and noted that 
transverse vibrations of a string also satisfy the wave equation and the same 
boundary conditions. 


3 The method of separation of variables 


You will recall that I have set myself the task of finding a function which 
determines the position of a vibrating spring at any time f. In other words, | have 
to find a specific solution of the wave equation 

au _1eu 

Ox* 


satisfying the boundary conditions 
U(0,1) = U(L, 1) =0 for all 1. 


It is clear that this equation as it stands does not contain sufficient information, 
for at no stage in the previous section have | introduced the initial position and 
velocity of the spring. However, as we shall see, if the initial position and velocity 
are specified we can find a definite solution to the problem. 


Before we get deeply involved in the details of the methods which we use to solve 
partial differential equations, | should like to say a little about the general strategy 
that we adopt. The point is that in general it is really rather difficult to find all the 
solutions of a partial differential equation. If, however, we restrict our search for 
solutions to those of a particularly simple form, then we may be able to reduce the 
problem from a difficult partial differential equation to a number of much easier 
ordinary differential equations. Such a process is the method which is known as 
separation of variables, and I shall explain the details of it shortly. 
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Let us suppose for a moment that we do have a method which will enable us to 
find some of the solutions of the wave equation 


@u _1#u 

a? oP” 
say U,, U2,.,.,U,; then since this partial differential equation is homogeneous 
linear, we know that any linear combination 


n 
YAU, 
fest 
(where A;, A2,...,A, are constants) is also a solution, from Theorem 1 on page 7. 
The method of separation of variables does indeed provide us with such a set of 
solutions, and by taking linear combinations we can obtain many more. But this is 
an idea which we shall develop later; at the moment our concern is to find some 
solutions of the wave equation. 


3.1 Separation of variables 
Our immediate plan is to find some solutions of the wave equation 
@u _1@u + 
a ar 
satisfying the boundary conditions 
U(0,1) =0 
U(L,1) =0 
We now fake ‘a leap in the dark’, guided only by the fact that other 
mathematicians, wiser than ourselves, have seen the final outcome, and attempt to 
find some solutions of Equation (1) by assuming that they take the form 
U(x, 1) = X(x)T (1). (2) 


In other words. U is the product of two functions X and T. where X is a function 
of x alone and T is a function of r alone. (There is of course at this stage no 
guarantee whatever that there are any solutions which take this form. By 
constructing them we shall in fact show that they exist.) 


hoe allt. 


Differentiating Equation (2) partially with respect to x, we obtain 
Ul) _ YeTt) 
Ox 


where X’ denotes the derived function of X. (Remember that 7, being a function of 
tonly, can be regarded as a constant for this calculation.) Differentiating partially 
with respect to x again, and abbreviating the notation, 


@U 

== =X'T. 

at x (3) 
Similarly, differentiating the original Equation (2) twice with respect to 1, we obtain 

@u * 

qe (4) 


Substituting Equations (3) and (4) into (1), we get 
pet yan 

X"T= qrAT “ 
Now dividing both sides by XT to get all the functions of t on one side and 
functions of x on the other, we obtain 

Be 8 

on 9% 
(assuming that neither X nor T is zero). This then is a necessary condition that 
U = XT should be a solution of the original equation. 


Do not confuse this use of T, 
to indicate a function of time, 
with our previous use of this 


letter to 
a spring. 


note the tension in 
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The next step caused me considerable difficulty as a student, and I needed to think 
about it for quite a while. We let 


ao es 
ae oT 
or 
x" an! 
bey (a function of x only), it must be independent of 1. But if y is independent 


then since p = (a function of ¢ only), it must be independent of x; but since 


of x and t, it can only be a constant. 


[I hope that I have made the point more clearly than it was made to me, but, if 
THX") 
OT) X(x)" 


you still find it difficult, let u(x,1) = 


Hx, = Lisl} = 10,1) (since each expression is independent of x) 
a T(t) 

(0,0) rey = (0,0) (since each expression is independent of r). 
Hence 


Hx, t) = 4(0,0) for all x and t, 
and therefore j is constant. } 


It follows that any solution of the original partial differential equation which is of 
the required form must satisfy the two ordinary differential equations 


T"=c'yT and X"=pX. (5) 


We know that j is a real (rather than a complex) constant (because T and T” are 
real), and therefore there are three cases to consider: 


u>0, w=0 and p<. 


Before I discuss the three cases in general, try the following exercises, the first to 
prepare you for the general argument and the second because the result will be 
useful later. 


Exercise 1 
Solve the equations 
T’ =cpT and X" =X 
in each of the following cases: 
@ c=2 p=l 


(i) c= w=-l 
(iil) e= -2,n=1. 
Exercise 2 


If h #0, and Ae“ + Be~™ = 0 for all values of t, show that A = B=0. 


[Solutions to Exercises I and 2 on p. 36) 


Remember that we are trying to solve the partial differential equation 

#u seu 

ox? 8 OF 
with U(0,1) = U(L, 1) = 0, and we have seen that a solution of the form 
U(x,t) = X(x)T(t) must satisfy T" = c? Tand X" = 4X. Now we use the 
boundary conditions to filter out the solutions which are of no interest in the three 
cases w > 0, w= 0 and p <0. 


The process is rather like panning for gold — we need to filter out quite a lot of 
useless material, until finally we are left with a nugget at the bottom of the bowl. 
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Case I: p > 0 

For convenience we can put 4 = h? where h is a real constant; then the equations 
T’=ch*T and X* =X 

have solutions 
T=Ae™+Be™ and X =Ce+De™ 

where A, B, C and D are constants. 


Now if either T(t) = 0 for all t or X(x) = 0 for all x, then U(x, r) = 0 for all x, t, 
and although this is in fact a solution of the partial differential equation, it is a 
rather boring solution since it corresponds to the spring remaining stationary. The 
solution U(x, t) = 0 is called the trivial solution. We are actually looking for non- 
trivial solutions which satisfy the boundary conditions 

U(0,1) = U(L,t) =0 for all 1. 


These boundary conditions refer to particular values of x (namely 0 and L) and all 
values of t, and assert that 


(Ce° + De®) T(t)=0 fopaale 
(Ce + De~™)T(1) = 0 a 


As we are searching for non-trivial solutions, we may assume that 7(1) # 0 for 
some value of ¢, and this allows us to conclude that 


Cc +D =0° 
ce + De™™ = 0. 


The Unique solution to this pair of equations is C = D = 0, which leads to the 
trivial solution after all, in spite of all our efforts to avoid it! 


Case I: p =0 
In this case the equations for T and X become 
T’=0 and x" =0. 


As in Case 1, all that really matters about T is our assumption that T(r) # 0 for 
some f. The solution of X" = 0 is 


X =Cx+D, 
and if we again apply the boundary conditions we obtain 
DT(1)=0 


(CL + D)T(t) = of forall’. 


Our assumption that T(r) # 0 for some f implies that 
D=0 
CL+D=0 
so that C = D = 0, and again the only permissible solution is the trivial solution. 


Case IL: xp <0 
For convenience we put 4» = —h?; then 

X" = ix and T= =c"T. (6) 
As in part (ii) of Exercise | we see that the solution of the first of these equations 
Is 


X = Ccoshx + Dsinhx. 

We now apply the boundary conditions (again simply assuming that T is not the 
zero function) to obtain 

CT( =0 


(CcoshL. + Dsin hL) T(t) = 4. Sa 
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and then choose a value of t for which T(t) # 0. We deduce that C = 0 and 
therefore 


DsinhL = 0. 


We may assume that D # 0 (for otherwise we are once again left with the trivial 
solution), and therefore 


sinhL = 0, 


This is the ‘gold nugget’ for which we have been searching. If sinhL = 0, then it 
follows that 


hL= nn, ie.h="™, 
for some integer n. 
Solving the second equation of (6), we obtain 
T= Acoscht + Bsin cht, 
which gives us the solution 
U(x,t) = X(x)T(D) 
= Dsinhx(A cos cht + Bsin cht) 


nm 
where h = 7 for some integer n. 


Substituting this expression for h and absorbing the arbitrary constant D into the 
arbitrary constants A and B, we finally have a solution 

ennt _ cnn 
Acos 5 ie Bsin er 
which satisfies the boundary conditions and which is certainly not the trivial 
solution. More than that, we have discovered an infinite number of distinct 
solutions corresponding to the choices of the integer n: 


MTX 
U(x, = sin | 


. AX cnt . CHL 
U(x = so (41008 = + Bsn") 

. 2nx 2ent . ent 
Ux) = o- (4,c0s 2 + Bz sin—— 


etc,, and the general function of this type would be 


mt 
——+ 3, si —— 


ne . nent 
L Bs 


U,(x, 2) = sin = [sn00s 


where n is an integer, and the As and Bs are arbitrary constants. 


These solutions resemble the normal modes of a lumped-parameter vibrating 
system, and what we shall do next corresponds to taking linear combinations of 
normal mode solutions to get the general solution. Although we have indeed found 
rather a lot of solutions of the wave equation, we have by no means found them 
all, and it will need some further work to increase our store of solutions. 


Exercise 3 
Use the method of separation of variables to show that if 


Ax, t) = X(x)T(t) 
is a solution of the partial differential equation 


then 
X*=cwX and Th =pT 


for some constant 


Unit 24 
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- (This partial differential equation may be used to model the flow of heat along a thin metal 
rod, in which case 0(x,1) represents the temperature at time / al a distance x from one end 
of the rod.) 


[Solution on p. 36] 


3.2 Initial conditions 

In the last subsection we discovered an infinite number of solutions of the wave 
equation (each satisfying the appropriate boundary conditions): 

nnet nnet 


—— + B,sin—— 


. AMX 
U,(x, t) = sin—— | A, cos L L 


3 
where n= 1, 2, 3,.... 
The current state of our investigation is summarized in the set diagram, where: 


= {all solutions of the wave equation}, 


# = {all solutions satisfying the boundary conditions U(0,1) = U(L,t) = 0}, 


= {all solutions of the product form U(x,t) = X(x)T(s)}- 


We have discovered that ¥ 7 #, the intersection of the sets ¥ and &, consists of 
solutions {U,(x, t)}; but our objective was to discover all the elements of the set # 
(shaded in Figure 1). 


We now use Theorem | (page 7), which tells us that we may add two solutions of 
the wave equation and obtain a third solution, but there is still the question of the 
boundary conditions. In fact if we add two solutions which satisfy the conditions 


U(0,t) = U(L,t) = 0 for all t, 


then we obtain a third solution with this property. For this reason these particular 
boundary conditions are called homogeneous. (We shall see an example of non- 
homogeneous boundary conditions in Section 4.) 


In fact any function of the form 
” 
U(x) = Y Unlx,0, 
net 


for some integer N, is an element of the set #. 


We are now ready to introduce the final factor into our investigation: 


The initial conditions 

In order to specify the motion of the spring, we need to know where it is and how 
it is moving at the start of the motion, say at time ¢ = 0. To know where it is we 
need the function U(x,0), which tells us the displacement of each point of the 


spring when t = 0, To know how it is moving we need the function BO) 
which tells us the velocity of each point of the spring when 1 = 0. 

au 
ot 
us to pick out the function U(x,1) from all the other functions in the set #, and 
thus to determine the motion of the spring. Let us now see how this works in 
practice, 

Example 1 


A uniform heavy spring stretched horizontally between the points x = 0 and x = 
is given an initial displacement 


If we are given U(x,0) and —(x,0), then this should be sufficient information for 


U(x, 0) = sinx 


along the length of the spring and released from rest. Find its subsequent 
displacement as a function of x and 1, if the wave speed (i.e. the constant 
appearing in the wave equation) is c. 


Figure 1. Various sets of 
solutions of the wave 
equation. 


The connection between the 
two senses of the word 
“homogencous’ (as applied to 
partial differential equations 
and as applied to boundary 
conditions) is that they both 
imply the property that ‘the 
sum of two solutions is 
another solution’, 
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Solution 

We can apply the wave equation 
@u _1@uU 
Ox? ar? 


for some constant c, and we try a solution of the form 


N 
U(x.) = YU, (x0) 
a=1 


N 
= ¥ sinnx(A, cos nct + B, sin nct) (7) 


Pos) 


for some integer N, where A, and B, are constants which have yet to be 
determined, 


We are told that 
U(x, 0) = sin x, (8) 


and since the spring starts from rest we have 
aU 
Fy 00) = 0. (9) 


Putting ¢ = 0 in (7), and using (8), we obtain 
sinx = A, sinx + Asin 2x + Aysin3x + «+: + Aysin Nx, (10) 


Then differentiating (7) partially with respect to 1, putting 1 = 0 and then using (9), 
we obtain 


0 = c(B, sinx + 2B; sin 2x + 3B, sin 3x + --- + NBysin Nx), (11) 


Our task is now to choose the As and Bs in order to satisfy Equations (10) and 
(11), and I have chosen this example deliberately so that it is possible to see these 
values at once. All we have to do is to choose 


A, =1, A, = 0 for n > 1, to satisfy Equation (10), 
and 

B, = 0 for all n, to satisfy Equation (11). 
Equation (7) now reduces to 

U(x, t) = sin x cos ct 
and this is the solution to our problem. 


At different instants of time t, the function U(x, 1) tells us the displacement of the 
point on the spring whose reference position is at a distance x from A. 


If we choose a particular point P on the spring, so that x is regarded as constant, 
then P moves backward and forward under the influence of the factor cos ct. The 

amplitude of the vibration is determined by the factor sin x, so that the movement 
is, as we should expect, most violent at the centre of the spring. 


Figure 2 (overleaf) shows the positions of various points of the springs at various 
times; for example P,, P;,...,P; represent various positions of the moving point P 
at times [),f2,....f7. 


Note the simplification that 
results in this example from 
choosing the length L to be 
equal to x, 
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— reference state 
ae a 


~d- . 4. 
+ + tn 
i 

+ + 4t0 


Figure 2. The movement of a longitudinally vibrating spring. 


S initial displacement 
4 
i 


It can be seen from this example that our task is finally reduced to finding the 
coefficients in two expressions similar to Fourier sine series, and thus at last we see 
the link between this unit and the preceding one. Using the techniques developed 


in Unit 31, we should be able to solve our problem in general. 


3.3 Finding the coefficients in the Fourier series 


Suppose that and y are given functions, each defined between 0 and L, with the 


properties that 
(0) = P(L) = WO) = Y(L) = 0, 
and we wish to find a solution of the partial differential equation 


@u 10U 
me 


on the interval 0 < x < L, which satisfies the boundary conditions 
U(0,1) = U(L,t) =0 for all ¢ 
and the initial conditions 
U(x, 0) = @(x) WO<x<L, 
au 


at (x0)=ix)  if0<x<L. 


(12) 


(13) 


Using the results of Subsections 3.1 and 3.2, we can attempt to find a solution in 


the form of an infinite series 


U(x, = = U,(x,t) 
, wet 


From Equation (12), and putting r = 0 in the above series, we obtain 


(x) = U(x,0) = ¥ A,sin 


ua 


L 


(14) 


(15) 
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From Equation (13), and differentiating the series in Equation (14) partially with 
respect to 1, we have 


ou S mmx [ onmc . . net | nme nact 
or De), an aaa A, sin L + L B, cos ate 
and therefore 
ou a ec 
Wx) = (0) = >} “Basin (16) 


nat 
We now use the techniques of Unit 31, Section 4, to determine the constants A, 
and B,, and to do this we use the functions and (which are non-periodic) to 
construct two new functions, ® and ¥, which (like the sine function) are odd 
functions and periodic. We let 

(x) = P(x) if O<x<L 


—o(—x) if-—L<x<0 


and 
Vix) = W(x) if O<x<L 
-W(-x) if-L<x<0 


and further require that both functions ® and are periodic with period 2L (see 
Figure 3). 


y= Mayor x) 
a 


Ee 
+ 
Pu 


Figure 3, Obtaining an odd function of period 2L from a function defined between 0 and L. 


We now use the formula for A, given in Subsection 4.1 of Unit 31, allowing for the 
fact that the period of our function is 2L, not L. We obtain from Equation (15): 


My ied _ NX 
An= [| Ooosin ax 
1/f nx be NRX 
= (/ O(x)sin a dx — f @(x)sin ma] 
alk )sit aes fou in™ as] 
=| neieteaeaset x)sin ds 


Lite MX a mnt 
= if (x)sin Tox + f $(c)sin + — de 
(putting r = —x in the second integral). 
Thus 
2 MX 
A= ral @(x)sin ora dx. 


Similarly we obtain from Equation (16): 


nme 2 rt _ MX 
7 B.=> |, (x)sin > —dx, 
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2. . Max 
B, -a), ¥(x)sin—— dx, 


To summarize, the general procedure for solving the wave equation with 
homogeneous boundary conditions is given in the box below. 


Procedure 3.3 
Given the wave equation 
au 1eu 
ox? a e ar 
subject to the boundary conditions 
U(O,t) = U(L,1) =0 for all t 
and the initial conditions 
U(x, 0) = p(x) fO<x<L, 


(O<x<L) 


oO s.0) = Hen) ifO0<x<L. 


1, Use the method of separation of variables to obtain solutions 


< WNX nnct ant 
U, (x, 0) = sin>— (4400s + B,sin ari 
2. Use the facts that the equation is homogeneous linear and the 
boundary conditions are homogeneous, to construct a general 


solution 
._ nnct 
L 


formulas for the As and Bs as follows: 


2 f° . MmX 
A= an (x) sin [—dx 


2 TX 
B= eI, Wo) sin™ de, 


« 
UxKo= ¥ sin (4400s + B,sin——}. 
Pan L 


3. Use the procedure in Subsection 4.1 of Unit 31 to obtain 


If the functions # and y are sufficiently simple, we may attempt to find formulas 
for A, and B, in terms of n; otherwise we may compute say the first twenty 
coefficients, perhaps using Simpson's rule and a modern computer to calculate the 


integrals. In either event our problem is essentially solved. 


Let us try this technique on a particular problem. 


Example 2 

Solve the partial differential equation 
@uU _1@uU 
ax? cat? 


on the interval 0 < x < x subject to the boundary conditions 
U(0,t)= U(x, =0 for allt, 

and the initial conditions 
U(x,0) = x(x — x) ifO<x<z, 


AO (ss0)=0 ifO0<x<x. 
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Solution 
We attempt to find a series solution 


U(x.) = ¥ sinnx (A, cosnct + B, sin net), Just as in Example 1, the 
eat assumption that L = x 
simplifies the formula, 


Applying the formulas in the above box for A, and B,, we obtain 
2¢" 
A, = =f x(n — x)sin nx dx. 
To 
Using the method of integration by parts, we can show that 
" 
f xsinnxdx = —"cos nn 
lo n 
and 
ere LZ | 
x? sin nx dx = ——cosnn + — (cos nn — 1). 
lo n n 
Thus 


2 2 2 
A, = 2{ Ecos + © cosnn - 2 (oanx - o} 
nl on n n 
4 
= wll — cos nz) 


4 
ert —(—1)} 


Since the initial conditions give y/(x) as the zero function, we have 
B,=0 for all n, 


and thus we have found a solution to the problem in the form of an infinite series: 


~ 


U(x, t) -4 yt — (=1)")sin nx cos net. 


This expression can be simplified, by noticing that 1 — (—1)" is equal to 0 for 
even n and to 2 for odd n, so that we have to sum only over odd n. Thus, putting 
n= 2m — 1, we finally obtain 


uxt =8 sin (2m — 1)x cos(2m — I)ct. 
(x0) =, Qm=1) a Pp )x cos ( Ye 
Exercise 4 
Solve the partial differential equation 

eu eu 

ox ce? or 


with the boundary conditions 
U(0,t) = U(x, t) = 0 for all ¢ 
and initial conditions 
U(x,0) = (x) ifO0<x<x 
aU (x0) =0 fO<xen 
et 
where (with d a given positive constant) 
2dx 5 n 
j= £0 Le 
d(x) = fO<x< 2 


=) er 
-u(i-) if><x<n 

x] 2 Figure 4. The graph of $(x). 
Hint: the solution to Problem 3(i) in Unit 31, Section 5 shows that for any positive /. 


2 __ nx 8d on 
7 [ otaisin ax =F 


26 MST204 323 


where 


(x) = 
[Solution on p. 36] 


Summary of Section 3 


In this section we have used the method of separation of variables to obtain a 
solution of the wave equation subject to certain boundary conditions and initial 
conditions. 

The method is summarized in the following flow chart. 
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Boundary conditions Initial conditions 


UW. = Lay = 0 


Wave equation 


U(x.) = dx) 
a 

Wem) = Wn) 
a 


for all, 


Look for solutions of the form 


Ul) = XO TU) 


Separation of variables gives 
T= CuTand.v" = aX 


for some constant 


4 
Low Ula) = (Cx + DITO) 
T= Acoscht + Bsineht pele where T (1) remains 
to be found 


X = Cooshx + Dsinhx 


Ulin = (Ce + DE) TU) 
where T (1) remains to be found 


The boundary conditions imply 
Ce Oand Dsink L =0, 
giving non-trivial solutions if 
hi = na/L. for some positive integer n 


‘Only the trivial 
solution obeys 
the boundary 
conditions in 
these cases. 


There are solutions: 


ieymua oF eter cr) 


‘The general solution with these boundary conditions is 


ute) = 5 sin F(A. cos + B, sin | 


I 


r 
x) sin = dx 


Use the fact 
that we can 
add solutions, 
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4 Transverse vibrations of a stretched 
string (Television Section) 


The following section may be read either before or after the television programme 
“Waves — Partial differential equations’. 
Suppose that an elastic string is stretched tightly between two fixed points A and 
B, and disturbed slightly in a direction perpendicular to the string; then the string 
will begin to vibrate. If U(x, t) denotes the displacement at time ¢ of a point P 
which was originally at a distance x from A, then U satisfies the wave equation 
#u _1 eu 
ax? oF 
for some constant ¢. 
The methods which we have developed for solving this equation can be applied 


with no change at all, and it is only in the interpretation of the solution that the 
difference from the case of longitudinal vibrations occurs, 


Example 1 


Suppose that the points A and B are a distance L apart, that U satisfies the wave 
equation, and that the initial conditions are 


U(x,0) = sin f0<x<L, 


0 (x0) =0 if0<x<L, 


Then using the same procedure as in Example | of Subsection 3.2 (page 20), we 
obtain 
a» RX net 

U(x, = sin -cos——. 
Figure 1 shows the position of the string at various instants of time, where 
P,,P2,...,P. represent the positions of the moving point P at times 1,,19,...,l6- 
Essentially then, this problem of the transverse vibrations of a stretched string is 
completely solved. 


New boundary conditions 

Let us now consider the effect of changing the assumptions in the problem, and 
instead of fixing B we shall assume that some external influence causes it to 
vibrate at right angles to the string in such a way that its displacement at time r is 
Rsinat. 

The function U still satisfies the wave equation, but we have altered the boundary 
conditions to 


U(0,1) =0 
forall. 


U(L,t) = Rsinat 


This is an example of non-homogeneous boundary conditions, for if we add two 
solutions of the wave equation which satisfy these conditions then we do not 
obtain another solution with this property. You should work through the next 
exercise to discover how this difficulty may be overcome. 


Exercise 1 

Assume that U(x,t) = X(x)7(t) is a solution of the wave equation, so that 
T' =CpT 
X" =x 


for some constant 


TV 32 


\\ 
Py 
\ Ps 
Pr 
Figure 1. The positions of a 


vibrating string at six instants 
of time. 
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(i) Use the boundary condition 
U(L,1) = Rsinot 
to show that 


Rsinwt 
TW) = Yap 


(ii) Substitute the above expression for T(t) into the equation T” = c?4T and show that 
2 
oy 


(iii) Solve the equation 
X"=yux 
2 
F oF 
with «=~ 
(iv) Use the condition U(0,t) = 0 to simplify the solution obtained in part (iii). 


(v) Write down a solution of the wave equation which satisfies the above (non- 
homogeneous) boundary conditions. 


[Solution on p. 36] 


At this stage you may well ask how we might find the many other solutions of the 
wave equation which also satisfy these boundary conditions. In fact we use a 
rather nice trick, and suppose that Uo(x, 1) is the solution which we have found in 
Exercise 1, and U(x,t) is some other solution satisfying these boundary conditions. 
Then consider the function 


V(x, t) = U(x, 1) — Uolx, 0). 


The essential point to notice is that Up and U satisfy the same condition at x = L, 
namely 


U(L, 1) = Uo{L, t) = Rsinwt for all ¢, 
and therefore 
V(L,t) = U(L,t) — Uo(L,1)=0 for all . 


Thus V(x, ¢) is just a solution of precisely the problem which we solved in 
Section 3, ice. it satisfies the wave equation and the boundary conditions 


V(0,t) = ViL,t) =0 for all ¢. 

Thus we can find V(x, f) as a series and then write 
U(x, 1) = Uo(x,t) + Vix, 

to solve our present problem. 


An example should make the method clearer, and for simplicity I shall choose 
L=n. 


Example 2 

Solve the wave equation 
eu _12U 
ox? 2 oP 


with boundary conditions 
U(0,t) =0 for all ¢ 
U(x, 0) = Rsinwt for all r 
subject to the initial conditions 
U(x,0)=0 for0<x<x, 


© (x0) =0 for0<x<z. 
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Solution 
We begin by letting 


OI 
Rsin— sin wt 

c 
Uk) =————> 
sin — 
ct 


which we know from Exercise | satisfies the wave equation and the boundary 
conditions but not the given initial conditions. 


Then we put 
V(x, t) = U(x, t) — Uolx, 0. 
Putting ¢ = 0, and using the notation of Subsection 3.3, we have 
V(x,0) = p(x) 
where 
(x) = U(x,0) — Uo(x,0) 
=0 forO<x <7. 
Also, 
M0) = Wx) 
where 
Wx) = 20 (x0) = 044.0) 


Rosin = 
=0= E 
. On 
sin —_ 
@ 
Using the box on page 24, we have 
V(x,t)= ¥ sinnx(4, cos net + B, sin nct) 
net 
where 
A, =0- for all n (as @ is the zero function) 


and 


B, = =f w(x) sin nx dx 
nace Jo 


—2Ro [*.. wx. 
= ——— _] sin—sinnxdx 
. On Jo o 


nne sin — 
€ 
os ; * 
sin{——n)}x sin |—+n]}x 
—R c . G 
? @ o 
ne sin — —-n +n 
c lo 


(using the formula for integrating sin xx sin Bx, on page 16 of the Handbook) 


Ro (2 +») sin(2— nx — (2 -n) sin (2 +1} ‘ 


a . OK lo lo 
ne sin — — +n} |——n 
c a ¢ ° 
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—- fo coerce de Sain 


nme sin 
‘ (F aie 
o 


2 cos sinnx + 2nsin = cos mx 


—Ro 


OT F J-= 
ne sin — 
c 0 
(using the trigonometric formula for sin (« + B) on page 8 of the Handbook, letting 
o 
a= == tn) 


—2Rwe _ On 
=——_——————_|sin— cos mn 
5 5\ cs ¢ 
n(«? — n*c*) sin 


on 
a7 Sin nx sin net 
nc? 


2Roc & (-1)"*! 
+ =——77 sin nx sin net. 
on Rye O* — re 


This then is the solution of our problem (which is referred to in the television 
programme), but you should notice that it is certainly not valid if sin =0,in 


other words if w is some multiple of c. In fact if @ is close to some multiple of c 
then the function 


«WX 
Rsin—sinwt 
Volx,) = 
_ On 
sin — 
i 


will be very large in modulus for some values of t. One of the terms in the series for 
V(x, 1) will also be very large in modulus, and near ¢ = 0 these two large terms will 
effectively cancel out, but as we progress in time, these large terms will dominate, 
and we should expect some very wild oscillations of the string. 


(There is a further difficulty which might have occurred to you. The point for 
which x = x (attached to the vibrator) is at a distance R sin wt from its rest 
position at time 1. Its velocity is therefore Rwcos wt at time t, and hence when 

t = it has a velocity Rew perpendicular to the string. But all of the adjacent 
points of the string are assumed to have zero velocity. There is no mathematical 
difficulty here, but there is certainly a physical one. To overcome this problem we 
might assume that there is a very small link, between the vibrator and the string, 
which is allowed to move, and whose effect we ignore.) 


The strategy which we have adopted to solve this problem and the physical 
interpretation of the solution are discussed in the television programme which 
accompanies this unit, ‘Waves—Partial differential equations’. 

Summary of Section 4 


The flow chart overleaf is intended to provide a summary of this section. 
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Wave equation 
®u_1 wu 
a 


Initial conditions: 
U0.) =0 


au 
—in0=0 
a 
INHOMOGENEOUS. 
BOUNDARY 
CONDITION 
U(L.t) = Rosin ot 


Look for a solution of the fotm 
Unde) = Xx) TH). 


Rein Separation of variables gives T’ = pT 
TO = —Soy— | | and. x” = wx (wconstant). 


Bring together 
the expressions for 
X(x) and T(t): 

Deancels out. 


Rain sino 
Uda) = 


‘The general solution with these boundary conditions is 
Ulx.t) = Up) + Vixen) 
where V(x,1) is the general solution of the problem 
with homogeneous boundary conditions. 


Use the flow chart on page 27 
to find V(x). 


‘The general solution with these boundary conditions is 
Rsin“ sine 
Ux) = =r + Disa 
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5 End of unit problems 


In this set of problems you will be led through an investigation of the vibration of 
a uniform heavy spring suspended vertically from a fixed point with a weight of 
mass M attached to its lower end. Later for simplicity we assume that M = 0. 

We use the notation of Example | of Section 2 (page 10), and suppose that the 
spring is of unstretched length /), modulus k and mass Mo, and that a mass M is 
attached to its lower end. 


The element between P and Q is presumed to have moved to become an element 
of length 6s at time 1, and P is then at a distance s from A. 
You may attempt this set of problems in the following form, or if you find it too 
difficult you may work through the sequence of questions given after the problems, 
which take you through the working in easy stages. 
Problem 1 
With the above notation and assumptions, show that 
@s | Mog _ Moé?s 

‘ae, a 
Problem 2 
Find a solution s(x,t) of the above partial differential equation which is independent of time 
(that is s(x, ¢) = X(x)). 
Problem 3 
Assuming now that M = 0: 
(i) show that the tension in the spring at B is zero; 
(ii) use Hooke’s law to show that 


és 
Allo = 1 for all t; 


(iii) find the general solution of the above partial differential equation in the form of a 
series, subject to the boundary conditions 
80,1) =0 } 
os for all 1. 
Allo) =1 


(Hint: Consider the difference between the required solution and the solution which is 
independent of time.) 


[Solutions to Problems 1-3 on pp. 36-8) 


Sequence of questions 
Our objective is to find a partial differential equation involving s(x, 1), the distance 
of P from A at time ¢, and we shall require the tension in the spring at P, which we 
denote by T. Read again through Example | of Subsection 2.1 (page 10) and then 
answer the following questions. 
(a) What is the mass of the element és? 
(b) What forces are acting on the element ds? 
(c) What is the total downward force acting on the element 6s? 
(d) Write down expressions for: 

(i) the distance of P from A, 

(ii) the velocity of P downwards, 

(iii) the acceleration of P downwards, 

at time 1. 


(e) Use the answer to Question (d) to write down an equation derived from Newton's 
second law 
force = mass x acceleration 


involving T and <>. Take the limit of this equation as dx —+0. 


@s 
er 


unstretched spring ~~ 


(with gravity B 
“switched off) 
stretched spring. 
(with gravity 
“switched on’) 


Figure 1. A spring in its 
unstretched and stretched 
states. 
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(f) Use Hooke’s law to write down an expression for T in terms of <. (See Subsection 2.1.) 
i 


. cf. es 
(g) | Use the answer to Question (f) to find an expression for x in terms of > 
é 


(h) Use the angwers to Questions (e) and (g) to write down a partial differential equation 
ai ee PH Re LL 
involving he and ra i.e. eliminate ax} 

(i) | What can you say about s(0,1)? 

(j) | What forces are acting on the mass M? 

(k) | What is the total downward force acting on the mass M? 


(I) What is the downward acceleration of the mass M (in terms of the function s)? 


22 
(m) Use Newton's second law to find a relationship between T(/o) and aill.t 


(n) What could you say about T(/o) if s were a function of x only? 


(0) What could you say about < when x =p. ifs were a function of x only? 


(p) Put s = X(x) in the partial differential equation obtained in Question (h), and hence 
integrate twice to find a solution which is independent of time. 

(q) Do you recognize the solution of Question (p)? 
(Look again at Example | of Section 2, page 10.) 

(rt) Make a substitution 

(M + Mog Mogx* 
k 2klo 


in the partial differential equation obtained in Question (h). What equation do you 
obtain? 


+1fx+ 


y= six) — 


(s) | What can you say about y(0.r)? 


{t) Assume now that M = 0. (The case M # 0 is too difficult for us at this stage.) What 
can you say about (lo), using the answer to Question (m)? 


(u) What can you say about Sit, using the answers to Questions (f) and (t)? 


(v) | What can you say about S(lo.t) using the formula for y in Question (r), and the 
answer to Question (u)? 


(w) Use the method of Subsection 3.1 (considering three cases exactly as before) to find a 
set of solutions of the wave equation obtained in Question (r) with boundary 
conditions 


0.1) =0 
for all 1, 


ey 
(bb. =0 
Fy lle) 


these being the boundary conditions derived as answers to Questions (s) and (v), 
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Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


14) F(xy)=— 


Pot 0b x 
x y wy 
(ii) F(x,y) = Ax? + 2Bxy + Cy, so 


oF 
ee 2Ax + 2By, 


(iii) F(x,y) = sin(2x + y), so 


oF = 2Bx + 2Cy. 
ay 


F 
ge = cos(2x + y). 


OF 
re 2cos(2x + y), oF 


(iv) Flxy) = 62x + y), 80 


F 
e = 29'2x +9) (where ’ denotes the 


derived function of #), 


"es (2x + y). 
(Y) F(x») = d(ax? + by?) so 


Gn = 2by p'(ax? + by*). 


ee ax (ax? + by), 
es a oy 


i) (8) 


i 
di) Sos0 
ax 


mn or 
(iii) eas (so 


al a 
3.41) F(x.y) = xy, 80 ue OF = x, and thus 
Ox ey 
OF oF 
oe Oyax ~ 


(ii) 


(ii) F(x,y) =e, so 
oF oF 


Be egy and thus 
OF @F @F OF _.., 
x? Brey OP 


Indeed, all the partial derivatives of F, of all orders, are equal 
toe”. 


4) U a fas so 


ou ue. 
e ¥ rise, = Fay 7PM 
eu eu 
US Ox Oy 


= Seg), 


Thus, = figs") 


eu 
oxdy" 
Gi) U =f(ay — bx), so 


ou 
ax 


= —bf"(ay— bx) and = aay ~ bx) 


eu ou 
Thus, at Uribe 


(ili) y = f(x — 1) + g(x + 0), 80 


® = poet gtx th 
ex 


ra 
= —six—tt eet) giving 
or 


ey 
Saale eH) and 


et =S"x— t+ ge + 0), 
ay ey 
=37 
5. (i) and (iv) are homogeneous linear. (ii) fails to be linear 


because of the term oe in which two partial derivatives 


are multiplied together, while (ili) fails to be homogeneous 
because of the term e*”, which does not involve the unknown 
or any of its derivatives. The only one of the four which is 
constant-coefficient is (iv). 


Solution to the exercise in Section 2 
1. Substituting s = U + x into Equation (13) gives 


e 


12 
me +x= aaa +3) 


Thus 
2U Ox 1 (Uo 
ax? “ax? ar © ot] 
ox @x ax ve) 

Now, == 1, 90 5 = 0. Also = =O since 5 means the 

tabesobscinsnge'os'? Varbes anelor reraae fied Thus” = 0, 


ar 
and we obtain Equation (14). 
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Solutions to the exercises in Section 3 
1,(/) Substituting c = 2, un = 1 gives 
T(t) =4T(t) and X"(x) = X(x). 
From Unit 6, Section 1, we have 
T= Ae™ + Be~™ for some constants A and B; 


also 

X =Ce*+De™* for some constants C and D. 
(ii) Substituting ¢ = 2, u = —1 gives 

T'=-4T and -X"=—-X. 
Hence 


T = Acos2t + Bsin2r 
and X =Ccosx + Dsinx 
for some constants A, B, C and D. 
(iii) Substituting c = —2, » = | gives 
T’=4T and X"=X, 


the same equations as in part (i). Thus the solution is the 
same as in part (i). This shows that although the sign of y is 
important, the sign of c is irrelevant. 


2. We know that Ae” + Be~™ = 0 for all values of t, so 
that, in particular, it is true for : = 0 and r = 1/h. Therefore 


A+B=0 
Ae + ae 0 
e 
and this pair of equations has the unique solution 
A=B=0. 
3. If 


(Xt) = X(x)T(H), 


then differentiating partially with respect to x (for fixed 1), we 
obtain 


0 
Ox? 


and differentiating partially with respect to ¢ (for fixed x), 


= X"(x)T(t) 


00 
a = X(x)T'(0). 


Tt follows that 
X"(x)T(t) = 7 X(x)T(0), 
so that 
Xe) _ TO) _ 
X(x) To) 


Using the same argument as before we deduce that y is 
constant, and hence the functions X and T satisfy 


X"=c*yX = and 


H, Say. 


T' =pT. 


4, Using the formula for 4,, we have 


Ay = f/ocesin nx dx, 
Rlo 
Using the hint (with / set equal to ), we obtain 
8d nn 


Ay = xsin=~ 
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Also B, = 0 for all n, since y(x) = 0, and it follows that 


in nx cos net 


ad 2 
URN =5 yr 
ot 

as (-1p"! 
Dy Dm =e 88 2m — I xc0s 2m — 1)er 


(This is the series used to model the motion of a vibrating 
string in the television programme ‘Waves — Partial 
differential equations’ which accompanies this text.) 


Solution to the exercise in Section 4 


1.4) Using the boundary condition 
U(L,t) = Rsint, 

we have 
X(L)T() = Rsint, so that 


Rsinwt 


Ti) = 7) 
Ma 77.5 


(1) 
(ii) Substituting the above expression for T(1) into the 
equation T"(t) = c?T, we have 

Ro*sinot _ uRsinot 


XL) XL) 
2 
so that p= > 


oF 


(iii) Putting p = — a 


in the equation X" = 1X, we have 


X(x) = Ceos* + Dsin. Q) 


(iv) If U(0,1) =0 for all ¢, then X(0) = 0, and using 
Equation (2) we deduce that C = 0, Thus we have 


XQ) = Dsin =. 3) 
©) Uet)=xX~TIN 
= Dsin®™ ee 


ol 
Dsin— 
c 


(using Equations 
(1) and (3)) 


Solutions to the problems and questions in 
Section 5 


The solutions to the problems are set out below as answers 
to the sequence of questions. The answers to Questions (a) to 
(h) constitute the solution to Problem 1; from (i) to (p) 
constitutes the solution to Problem 2; and from (q) to (w) 
forms the solution to Problem 3. 
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1.(a) The element ds was originally of length 5x, so that its 
Modx 


mass is 
lo 
(b)(i) The tension upwards at P; 
(ii) the tension downwards at-Q; 
(iii) the weight of the element ds acting downwards. 
Mogox 


(©) Tee + x) — Tx) + =F 
lo 


i) six). 

3. Os “ és 

(ii) Fra) {sbbreviaea to < 

Peas) é es 

(iii) ae) {abbreviated to 3). 

Mogdx _ Modx @s 
lo I, or” 

and in the limit, after dividing by 6x, 


eT | Mog _ Mo és 
Io or?” 


(e) T(x + dx) — T(x) + 


x 
6s 
ne a — Sf 
qf) T ax i 
(Use Equation (3) on page 10 with T, = 0, since the spring is 


originally unstretched; and use a partial derivative because s 
is a function of x and 1.) 


oT _ as 
@) aha 


@s Mog Mo és 
a Ak 


This is the solution to Problem 1. 


2.) The top of the spring is fixed, so that s(0,1) = 0 for 
allt. 


(j) Its weight Mg downwards and the tension in the spring 
upwards. 


(kK) Mg — Tilo). 
2, 2 
(} SF losts in other words oe evaluated at (I, 1). 
a 
(m) Mg ~T\lo)=M—(lo.t). 
eid as @s 
(n) T(lo) = Mg, since in this case Fd would be zero. 


és COPae ‘ 
(0) Gx Cam be replaced by — if's is a function of x only, 


and from the answers to Questions (f) and (n) we have 


Putting x = /g and using the answer to Question (0) we have 
Mg + k= —Mog + A, 


so that 

A=(M+ Moje +k 
and therefore 

AGE = (M + Moje +k MOE 
Integrating again, 

Kx = (M+ Mode + Wx — MOB 


and from the answer to Question (i) we have B = 0, so that 
_ [M+ Mole Mogx? 
x= (Maes i) kl” 


This is the solution to Problem 2. 


3.4q)_ This is precisely the solution which we obtained in 
Example | of Section 2 (see Equation (8) on page 12). In this 
context it is just what we should expect—namely the 
solution which is independent of time. 


(tr) The wave equation, 


This is the hint suggested at the end of Problem 3. 
(s)_ 9(0,1) = s(0,1) = 0 for all 1. 

() T(lo) =0 because M = 0. 

This is the solution to Problem 3(i). 


(u) From the answers to Questions (f) and (t), 
és 
Rod = 1 for all ¢. 
This is the solution to Problem 3(ii). 
(v) Using the formula in Question (r) (and putting M = 0), 
ay 2 _ (Moe , |) , Moss 
Gx 6x k kly * 
so that using the answer to Question (u), 
Mog Mos 
ic % i a: 
=0 forall. 


(w) The discussion of Cases I and I] is almost identical to 
that of Subsection 3.1, and there are no non-trivial solutions 
of either form. For Case III we have 

v(x, 1) = (Ccos hx + Dsinhx)T (0). 


Since y(0,t) = 0 for all t, we must have C = 0 for a non- 
trivial solution. 


oy 
3x lor) =1 


Since 2 (lo, t) = 0 for all 1, we have 


Dh cos hl, = 0, 
so that 
hl = 5 + Ld 


for some non-negative integer n. 
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The required set of solutions is therefore 


valet) = sinin +4)" | Aycosin + $5 + Bysin(n +4) 
fo a 


This is the solution to Problem 3(iii). 


The gerieral solution for y(x, 1) is therefore 
¥ sinn + 4) 
n=O lo 
and so, undoing the substitution of Question (r) while 


remembering that now M = 0, the general solution for the 
original function s(x, 1) is 


M, Mogx? 
soars (Mat +1) x= “a 
lo 


(4, cos + rRaaal + B,sin(n + Rig 
lo Iy 


+ ¥ sinon + y= (4,cos +t + By sin(n + 4)". 
=o Io bo lo 
As in Subsection 3.3, a knowledge of the initial conditions 
would allow us to calculate the values of the A,s and ,s, 
thus solving as a Fourier series the problem of a heavy 

spring vibrating vertically under the effect of gravity. 


